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Introduction 

Let □ denote a wave operator on M x M'^, d > 1 of the form 



where A(x, Vx) is a time- independent, uniformly elliptic partial differential operator of order 
2 on R'^. 

As well known, if v : M+ x M'^ — )• M solves the wave equation with initial condition 



(0.1) 



nv{t, x) = 



v{0,x) = 0, dtv{0,x) = F{x) 



(t>0) 



(0.2) 



where F : M — )• M'^ is a smooth function rapidly decaying at infinity, then, letting v{t,x) = 
{E,{t)F){x), 

w{t,x):= / {Ei{t- s)F{s,-)){x)ds (0.3) 
Jo 

solves the equation \Jw{t, x) = F{t, x) on M x W^. We write w = G{F) and call G the Green 
kernel of □. 

The problem we are interested in, and solve in the case of the Klein-Gordon operator on 
Minkowski space-time, is the following. Can one find a n-independent 'optimal' Riesz basis 
('0A)AeiD) of L^(M X M'^) [12j, indexed by some set D, in which the matrix coefficients of the 
Green kernel, 

Ga,A' := (V'A,GV'A'), A,A'gD (0.4) 

are essentially diagonal, with fast off-diagonal decay ? 

The principal motivation for this work comes from the rigorous study of real-time quan- 
tum field theory [8t 120 ^ \T9 \ \10 \ \TT\ . From the beginning, constructive methods in QFT have 
been applied to imaginary-time, i.e. Wick-rotated theories, thereby replacing hyperbolic 
wave operators by elliptic operators whose study is much simpler. Alas, real physics lives in 
real-time and is founded on perturbations of linear, hyperbolic operators: Schrodinger oper- 
ators describing quantum many-body Hamiltonian systems) , or wave- like operators describ- 
ing relativistic particles of high-energy physics [TTj. The relevant quantities, including the 
simplest of all, the Green kernel, are then oscillatory integrals which cannot be appropriately 
bounded pointwise but only in average. Thus this article gives a satisfactory representation 
of the Green operator for the Klein-Gordon operator on Minkowski space-time, a good start, 
hopefully, for the constructive study of quantum field theory on pseudo-Riemannian man- 
ifolds. We also hope that such bases eventually turn out to be useful in the mathematical 
study of nonlinear wave equations (global well-posedness, blow-up, global parametrices, see 
e.g. [23] or |26j for an old review) or of dispersive estimates for linear wave equations with 
variable coefficients, such as the wave equation on curved space-time, or on domains with a 
boundary j3l[l3l[l5j. 

A very similar problem has been studied by Emmanuel Candes, Laurent Demanet and 
David Donoho |4|15j a few years ago. They were looking for an optimal representation of the 
solution operator E{t) : mo(") ^ '^^{t-, ) of linear symmetric systems of first-order, hyperbolic 
differential equations of the form 

d d \ 

- + J^^fc(x)— + S(x)j u = 0, ^/(0,x) = uo(x) gR". (0.5) 

k=l / 

The Riesz basis they constructed is made up of generalized, highly anisotropic wavelets 
for which they coined a nice name, curvelets . The matrix representation of the solution 
operator is not diagonal, but rather similar to a permutation matrix, curvelets being shifted 
in phase space according to the underlying semi-classical Hamiltonian flow obtained by 
the geometric optics approximation. Candes et al. wrote a series of papers starting from 
2004 on the remarkable properties of these curvelets, both theoretical studies (representa- 
tion of Fourier integral operators [6] and applications to the mathematical representation 
of images with edges; curvelets are now commonly used for algorithms of compression of 
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numerical photographies, which are much more efficient than those based on the conven- 
tional, isotropic wavelets. The interested reader may refer to the web-page of E. Candes 
(|http^www-stat. Stanford. edu/~ candes/publications.html) for more information and refer- 
ences. Curvelets may be seen as a phase-space extension of the second Paley-Littlewood 
dyadic decomposition used by Fefferman [9], Seeger, Sogge and Stein [21] to obtain esti- 
mates for Fourier integral operators. More recent references are two papers by Hart F. Smith 
[22| [23] (followed by a series of papers by Tataru [271 EH EH] relying on a totally different 
method); in the first reference. Smith proved Strichartz estimates for wave operators with 
C^'^ coefficients by introducing a kind of curvelet decomposition which partially inspired 
the work by Candes et al. 

The main feature of curvelets is that they are directional, anisotropic objects with a 
parabolic scaling, namely, their 'size' is 0{2~^) in one direction, equally distributed on the 
unit sphere, and 0(2~-^/^) in the d — 1 transversal directions; 'size' refers to the bulk of the 
support of the curvelets since they are not compactly supported. Then matrix coefficients 
of the solution operator E{t) have a fast decay off the shifted diagonal in terms of an 
appropriate pseudo-distance between the bulks of the supports. 

Of course, scalar wave operators may be converted into 2 x 2-systems of first-order, 
hyperbolic differential equations, and our operator Ei{t) is essentially a particular case of 
E(t). However our problem differs, in that we are interested in the Green kernel, given 



by the time convolution by Ei as in eq. (0.3). Whereas the curvelet basis that Candes, 



Demanet and Donoho constructed are functions of space coordinates only, we define a basis 
made up of functions of space-time. Quite naturally - in view of our application to wave 
operators we called these new objects Minkowski curvelets. 

Minkowski curvelets have three scales instead of two, given in terms of two independent 
scale indices m,j called respectively main scale and secondary scale. As in the case of 
curvelets, they are not compactly supported, but they have a fast decrease outside a small 
space-time window corresponding to the 'bulk' of their support. The space-time windows 
have largest size rays (longitudinal direction), corresponding to quasi- 

classical propagation in phase space. Their size along the independent transversal, spatial 
directions is 0(2^*"). Finally, going away from the direction of propagation in an arbitrary 
off-propagation direction, curvelets fade away after a distance 0(2^™~-'/^). Restricting to 
spatial directions and rescaling by typical propagation distances 0(2"™+-'/^), one obtains 
the same ratios 0(2"-'), 0(2~-'/^) as for curvelets. 

Our main result is the following. We let (V'a) AeD be a fixed basis of Minkowski curvelets. 
To each index A are associated in particular two scale indices, the main scale rriA and the 
secondary scale j'a- We define in section 2 a scaled 'distance' d{A, A') between two curvelets 
with indices A, A'. 
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Main theorem. 

Let G he the Green kernel of the Klein-Gordon operator := — X]^=i ~ A*^ 
Minkowski space-time, with mass /i > 0. Let A e D™J,A' G 'J' , with m,m' G 
j, j' > 1, m + j/2, m' + j'/2 > 0. Then, for every N > 0, there exists a constant Cn > 
such that 

\{Gi^A,M\ < 2-2-P(-'-')rf(A, A')-^. (0.6) 

This result holds true with a slightly different definition of the Minkowski curvelets, adapted 
to the Klein- Gordon operator. 

Assuming = "t-a' = and j'a = Ja' = J ) and using the relativistic notation xq for 
t, Fourier variables jei)*=o,..., 3 = (Cc^) scale like ||eo| - 1^1 1 ~ 2—^/2, |eo| + |^| 2^+^/2 
on the support of V^AjV'A'; so — ^2 _ ^2 ^^j^^ symbol of Do) scales like 2^"^ and a 



rapid estimate yields G p « 2 , accounting for the scaling factor in eq. (|0.6|). The 



interesting point is that the matrix Ga,A' has an almost exponential off-diagonal decay. 
In other words, the Fourier cut-offs destroy the propagation of the waves, at least in the 
ultra-violet region m + j'/2,m' -\- j'/2 > where |^| is large. 

Our estimates in the theorem are optimal in a region l^oM^I ^ !> ^ ^ 1 where 
geometric optics is a valid approximation. They are not optimal (or even wrong) in the 
irrelevant regions |^o| ^ 1^1 ^ 1 or |^| ^> |,^o| ^ 1 corresponding to j = 0, and in the 
infra-red region m + j /2,m' + j'/2 < where |^o| and |^| are not large enough and terms 
of lower order in the wave operator become essential. A companion article in preparation 
is devoted to complete, optimal estimates of the Green kernel of the Klein-Gordon operator 

on M X W^. It requires the introduction of yet another family of curvelets, Schrddinger 
or parabolic curvelets, for the study of the infra-red region, corresponding to the bottom of 
the mass hyperboloid |Co| ~ Z^; 1^1 ^ A*- 

It is reasonable to hope that the same type of result (or maybe a weaker one, maybe with 
a restriction to a subset of scales) is true for more general wave operators, such as e.g. the 
Klein-Gordon operator df — A — ff^ on M x M, where AI ~ M'^ is equipped with a Riemannian 
metric which is asymptotically flat at infinity jl5j. Preliminary computations show that a 
naive short-time expansion of the solution map Ei{t) using Fourier integral operators does 
not give precise enough estimates, even for bounded times, if □ is a general wave operator. 
Obviously, the precise definition of the curvelets should depend on the operator itself, since 
this is already the case for the Klein-Gordon operator on flat space-time. Of course one 
expects that the basis of curvelets does not need to be changed for 'weak' perturbations of 
a given wave operator. The difficulty is that a semi-classical analysis is not sufficient to get 
precise estimates. 

Our construction relies heavily on previous work by Candes, Demanet and Donoho, in 
particular on reference The article is organized as follows. We define our new objects, 
Minkowski curvelets, in section 1, extending en passant the construction of usual curvelets to 
arbitrary dimension (originally they were introduced in dimension 2, essentially for simplicity 
of exposition). Properties of Minkowski curvelets, in particular precise decay estimates, are 
developed in section 2. Section 3 is devoted to the Main Theorem. 

Notations. If t G M, (t) := 1 + |t| and {{t)) := 1 + \t\ + \t\~^. (x-')jez is a smooth 
multi-scale Fourier partition of unity (see Appendix). Minkowski signature is by convention 
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(+,-,-,-), i.e. x-^ = xo^-xi^i-xa^-a^s^ andx^ = (x^)'^ -xj-xl-xj, = ^2_^2_ 
^1 — The spatial components of x, resp. ^ are denoted by x, resp. ^; the reader should be 
warned against mistaking x"^ (Minkowski squared norm) for |a;p (Euclidean squared norm of 
spatial components), while \x'^\ is simply the absolute value of x^. The Fourier transform is 
defined with respect to the Minkowski scalar product, namely, Ff{^) = /(^) = / e^^'^ f{x)dx. 
The symbol x < y means: \x\ < C\y\ for some universal constant C. The symbol x y 
means: C~^|x| < \y\ < C\x\ for some universal constant C > 1. 



1 Curvelets: a reminder 



Candes-Donoho-Demanet [U [5], relying on estimates for parametrices of wave equations 
previously developed by Hart F. Smith [22], introduced a family of generalized wavelets 
on which they called curvelets. As we mentioned in the Introduction, curvelets are 
appropriate for the study of parametrices of wave equations (i.e. for the solution of wave 
equations with initial conditions), in that the matrices e'*"^ := ((e^*^* "^^'•^^^^fe V'Aj '0A'))a,A' 
are maximally sparse, i.e. {e^^^)A,A' is maximal when A' = A(t) is the image of A by the 
underlying semi-classical Hamiltonian flow, and is fast decreasing in the distance between 
A' and A(t) otherwise. 



1.1 About wavelets and frames 

By wavelets we mean orthonormal wavelets. The profane reader may refer e.g. to \12\ [T5] 
for a simple and pedagogical introduction to real harmonic analysis and wavelets, and look 
up in the classical references [HI [25j for detailed proofs and more advanced stuff. Let us 
simply recall here the following definitions and properties of wavelets. 

A (one-dimensional) wavelet is a function t/j G L^(M) such that the functions {il^-''^)j^-L^k& 
defined by Fourier rescaling and scaled translations in real space, V^^''=(x) := 2J'/2^(2i(x - 
k2~^)), make up an orthonormal basis of L^(M). 

The following two equations characterize wavelets (see [12], chap. 7): 

Y,m-'o\' = h (1.1) 



^'0(2-^'^)V;(2-i(^ + 2m7r)) = 0, ^ G M, m G 2Z + 1. (1.2) 



We shall only use eq. (1.1) that we shall call the first basic wavelet equation. 



There are many different constructions of wavelets. The easiest ones are obtained by 
multiresolution analysis through Fourier analysis, and are compactly supported in Fourier 
space. Assume, say, that [suppV'l C [1/2,2]. Note that i^^'^'iC) = 2-^/'^e-''''^~'^{2-^ 
hence |supp 

^ [2i-\2^+^] independently of the space localization index k; then j plays 
the role of a momentum (or Fourier) scale. Lemarie and Meyer constructed a family of 
wavelets along these lines, depending on the choice of a more or less arbitrary 'bell function' 
satisfying some support and symmetry properties, and satisfying the above support property. 
Provided the bell function is smooth, ip is compactly supported and smooth, hence ip is fast 
decreasing at infinity, i.e. \iIj^'^{x)\ < Cr{x — k2^^)^^ for every r > 0, with a constant 
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Cr depending on r. In other words, is centered around a window (an interval in 
fact), say, [{k — 1)2"-', {k + 1)2"-'], of size ^ 2"-', and has fast decrease outside this window. 
Characterizing a wavelet {'ip^'^) by its window A-'''^ := [{k — 1)2"-' , [k + 1)2"-^] , one may index 
the orthonormal basis of wavelets by the windows, {ip^''^)j^k = (V'a)agDj where D = Wj^jJ])^ 
and W is the set of intervals of length 2^~^ centered at points which are integer multiples of 
. This is the type of wavelets we shall be using. Later on we shall introduce generalized 
wavelets {curvelets or Minkowski curvelets) for which the couple of indices {j,k) becomes 
some more complicated set of indices, and it will be more convenient to write the wavelet 
family as (V'A)AeD5 where A is a window which is some compact subset of M'^, and A ranges 
in some complicated set D. 

For many purposes, it is not really necessary that {ij)^'^) constitute an orthonormal 
basis. It is enough that they should make up a tight frame (or frame for short in the 
sequel). A frame {f^)j£i, I ~ Z, is a family of elements of a Hilbert space 7i such that 
the Parseval identity is satisfied, i.e. = J^ji^^ f^)'^ ^ equivalently (see [12], Theorem 

7.1.8) (f) = ^j{4>-, f-') f-' for any cj) £ Ti. Orthonormal bases are automatically frames, but the 
converse is false; the frame property means simply that the sum of projectors 'Ylji'-i f^) 
is the identity operator, which is possible even when the f^ are not linearly independent 
(typically, 1, e^''^/^, e^^'^l^ G make up a frame of M^). 

One checks immediately that the Fourier transformed family (/■')jg/ is also a frame if 
n = L2(M^). 

Higher-dimensional wavelets (say, living on M'^, d >2) depend on the choice of a scaling 
matrix A; if they are compactly supported on Fourier, their support must be an enlargement 
of one of the fundamental domains of M'^ under the action of A] in practice, say, |^| E 
[2-'"^, 2-'"'"^] ii A = 2Id. In this case there is only one momentum scale. On the other 
hand, in our case, we are interested in generalizations of wavelets having several different 
momentum scales. In the case of the construction by Candes-Donoho-Demanet, the norm 
scale j corresponds to |^|, and the angular scale j /2 to the direction or angular sector cut 
along the sphere |^| = 2-^. A direct tensorization of the one-dimensional construction is not 
possible since the angular scale depends on the norm scale. The authors obtain instead 
frames which are almost as good as actual wavelets, in the sense that the scalar product 
of two elements of the frame have a fast decrease in the 'distance' (to be defined) between 
their windows. 

As is often the case, we shall disregard the infra-red region j <0 and 'resum' the (■0"')j<o 
into a 'father wavelet' ■00 defined through its Fourier transform. 



V i>i 

1.2 The Candes-Donoho-Demanet curvelet construction 

The following lines are not a direct repetition of the contents of [H [5] since we are in- 
terested here in three space dimensions, but obviously the passage from 2d to 3d is very 
straightforward. The construction relies mainly on two fundamental remarks: 
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Lemma 1.1 (i) (scaled Fourier window property) Assume that (j > 1,— 2~-5'/^ < 

I < 2-'/^) is a set of L'^ -functions such that for all ^ G R^, = 1; ^-fi- 

satisfying a scaled Fourier window property. Then 'Yji |V'"'(?)^"''HOP = J2j IV'"' (OP- 
In partic ular, if ij) is an orthonormal wavelet, hence satisfies the first basic wavelet 



equation (1.1 



), then also satisfy the first basic wavelet equation. 



(a) (space localization trick) Assume (/-')jgN is a family of functions in L (M) such that 
supp(/-^) C [0,2-K ■ 2^] and, for all ^, |/-^(C)P = 1, i.e. satisfying the first basic 

wavelet equation. Then /•'''^(x) := 2~^l'^f\x — k2~^), j £ N,k £ Z make up a frame 
o/L2(M). 

Proof, (i) is trivial. One way to construct such functions P'' is to start from a smooth 
function i> : M — ^ M with support in [— 27r, 27r] such that 11^(01^ + \^{^ ~ 2vr)p = 1 for 
< ^ < 27r. Then an easy computation shows that Yl-23/^<i<23/2 I'^'^'HOP = 1 on [0,27r], 
where i>-?'^(0 := i>(2J/2(^ - 27r • 2-^/^1)) is supported on 27r/2~j/2 + [-27r2-j/2, 27r2~-'/2]. 

For (ii), note that (0, /■''''')l2(r) = {4>f,'^~^^^e-'''^~'^)Lmo,2n-2^]) by the support condi- 
tion; thus, {(f>, P'^)l'^{m.) is the k-ih. Fourier coefficient of (pp viewed as a 27r • 2-' -periodic 
function. Using the usual Parseval theorem for Fourier series and then the first basic wavelet 
equation, one obtains 

Ei^'^'/''')|' = E/ \^Pio?di = I \m\'d^- (1-4) 

j,k j 

Hence {P'^)j± make up a frame of L^(]R). 

Thus one easily cooks up a recipe for constructing a frame of "curvelets" (V'a) on as 
follows. Recall from the preceding paragraph that there are two momentum scales, the norm 
scale j, and the angular scale j'/2. The windows A are now rotated rectangles of the form 
aiA + 2"-', 2"-'] X [— 2"-'/^, 2~-?/^]^), where Vf^ is some rotation sending the vector 

(1, 0, 0) into the unit vector called direction of the curvelet. One requires that (i) the ca, 
A £ W are regularly scattered 0{2^) vectors on the unit sphere 5^ := £ M'^ | |,^| = 1}, 
namely, |eA — e^'l ~ 2~-^/^ if A' 7^ A; (ii) the Fourier support of the curvelets is 'dual' to 
their support in direct space, i.e. (7^^'^)-l(supp V^a) C [2^-^,2^+^] x [-2^/^ 2J/2]2_ jj^^^^ 
a possible construction of a frame (V'a) satisfying the required properties: 

(i) (radial dependence) Start from a Lemarie-Meyer wavelet on M, ip. 

(ii) (angular dependence) Introduce two local maps on the sphere 5^, Map: U ^ (— 7r,7r)^ 
and Map : tj — )■ (— 7r,7r)^ (left composition of Map by the symmetry (xi,X2,X3) 1— )• 
(xi,X2,— X3)) defined on a neighbourhood U, resp. U = {(xi,X2,— X3) | (xi,X2,X3) £ 
[/} of the North (resp. South) hemisphere in 5^, and two functions /, / : — )• 
M4- such that supp I C U, supp I C U and P + P = 1 . These maps will al- 
low us to define curvelets with direction centered around unit vectors e = ei-^^j^ = 
Map~^(7r/i2-j/2,7r/22-j/2)^ _2i/2 < /^^/^ < 2^/^ or the symmetric ones ei^^i^.' Let 



then i)-'''^'i''2 (^j^j ^2) = y-''^^ {9i)v^'^^ {62) be defined as in the proof of (i) in Lemma 1.1 
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and 

= V.(2-'|^|)/(|^)i>^-'-(Map(||)), resp. = V.(2-'|^|)/(|^)i>^-'-(Map( ||)) 

(1.5) 

By Lemma 1.1 (i), the functions (V'-''*^, V'"''*^)i,e satisfy the first basic wavelet equation. 



In the sequel, e, Map, U shall stand either for e = e;^ Map, U or their symmetric 
counterparts e, Map, U. 

(iii) (space localization) Fix j and e. Then ■0-^'^ o 7^*^ is supported in the anisotropic 
rectangle [0,2-'"''^] x [— C2-'/^, C2-^/^] for some large enough constant C. By a simple 
rescaling one may assume that C = 1. Applying the space localization trick (see 
(ii)) to the three space coordinates yields a frame (V'"''^''^)j,e,fc) G Z^, 



Lemma 
withQ 



1.1 



iP^'^'f^ix) = 2-H^'^{x - n''{ki2-^,k22-^/^, hi2-^l'^)). (1.6) 



By construction the Fourier supports of the curvelets essentially do not overlap: {ip^'^'^^ip^''^''^') 
except if j and j' , e and e' are Fourier neighbours, i.e. |j — j'| < 1 and |e — e'| < 2~^l'^ , a fi- 
nite, j-independent number of possibilities for (j, e) fixed. On the other hand, the amplitude 
|^i,e,fc(2;)|2 is 'maximal' on A = A{j,e,k) = n''{{ki2-i ,k22-i/'^ ,k^2-il'^) + [-2-^,2-1] x 
[—2^-'/^, 2"-'/^]^), a thin, flat piece of wall (for j large) perpendicular to the curvelet direction 



Again by construction, |V'a(^)| ^ 2 ■'^ |^Vol(supp V'a) j and (as a simple variable 

rescaling shows) \i^/^{x)\ < 2^^ ^ (Vol(A))^5. Note that Vol(supp ^a)VoI(A) ?s 1, an 
optimal upper bound in view of the uncertainty principle. 



1.3 Minkowski curvelets 

The previous curvelets are a nice tool for the study of the solution of a wave equation at fixed 
time t. On the other hand, controling e.g. ||G(-F)||^P(^g) in view of Strichartz estimates 
requires the introduction of a new type of space-time curvelets (^a) essentially supported 
inside a compact space-time window, and such that G{t,x;t' ,x') is very small when x and 
y belong to far away windows. This is a priori contradictory with the fact that waves 
propagate. In order that G{t,x;t' ,x') be small for x' ~ x[t) following the semi-classical 
Hamiltonian flow, one must impose a Fourier support condition on which contradicts 
propagation. 

Let us proceed a little informally at this stage. Replace □ by the usual wave operator 
— A, so that G is the convolution kernel 

G{t,x) = I '-^e'-M, (1.7) 
see e.g. [51], or (in relativistic coordinates) 

G{t,x)='-}^j 5{io-\$.\f-j^di- y 5(eo + iii)^de|. (i.s) 

^Note the scaling prefactor 2^-' = ^2^J x (2^^/2)2^ jjj ,^ two-dimensional setting [3] one would find 
instead 2-^^^^ = V2--' x 2-^/2. 
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The latter expression emphasizes the fact that the essential contribution to 
G(£) = lim^.n+ ^. — comes from the region fn = ±1^1- Assume for the sake of 
the discussion that > 0. In order to contradict propagation, it is natural to require that 
^0 — 1^1 be bounded away from 0. Thus one is led to introduce two scales m G Z and j G N 
and a corresponding smooth cut-off X™'"'(0 which 

leo - 1^1 1 ~ 2™-^'/2, ^ 2'"+J'/2. (1.9) 

In the following discussion we assume j > 1 so that |^o| ~ 1^1- These constraints are 
preserved by transverse variations (5^_|_ _L ^ of ^ of order 1 6^j_ \ < 2™ since 

1^ + - m = 1^1 (V^^WINW) « l^-^ < (1.10) 

We now consider the bulk of the sup por t of G in space-time. Assume to simplify that 
X is parallel to i. Then e"'^'? in eq. ([l8| may be rewritten as e"'[(^''"l^l^^°+l^l(''°"l^l)] . 
Oscillations due to the phase imply fast decay of x'"'"' * G* x™'''' when \'3^\ <C ('^o — |^|) ""^ ~ 
2-m+jr72 — ^ l^l""*^ ~ 2~™~-^/^. In particular, this means that waves propagate 

only up to a propagation time of order 0(2"™'^-'/-^). 

The construction of Minkowski curvelets obeys the following principles. Since 1(^^x1/1^1 ~ 
at most, it is clear that a proper Fourier space decomposition should include an 
angular decomposition of ^ with 0(2-') sectors, in precise analogy with the case of curvelets, 
up to a rescaling by the lower Fourier scale — 1^1 1 ~ 2™"-'/^, namely, 2^ and 

Ig^^l^ll ~ 2-'/^. Correspondingly, the shape of a Minkowski curvelet in real space should be 
similar to that of curvelets up to rescaling by \x\ i.e. the space-time coordinate dual to 
Co-|^|. 

Let us now be a little more precise. The inverse Fourier transform yields the exponential 
e~^^'^, with X • ^ = xoCo — x ■ Assume the direction of propagation is along the first 
axis of coordinates, ca = (1, 0, 0). Then x • ^ « xq^o — = 2;o(^o — ^i) — (^^i — xo)£,i ^ 
xoiCo - 1^1) - - xo)\^\. Since - 1^1 ~ 2"^-j/2 and ~ \C\ » 2""+^/^, one obtains 
by dualizing |xo|,|a3| < 2-'/^"™ and |xo — |a^|| < 2^-'/^"™'. Assume \x\ ^ 2-'/^"™', which is 
generically the case. The constraint |xo — l^cH ^ 2"-'/^"™' ^ \x\ ~ 2-'/^"'" imposes (as in 



eq. (1.10) above) \xi\ < 2~"^, i = 2,3 for transverse coordinates, dual to the condition 
I'^il ^ 2™. In other words, rescaling by the typical time of propagation |xo| 2-5'/^"™' 
and replacing xi by the proper invariant coordinate along light rays, xi — xq, i.e. letting 
7/1 := 2™'~-'/'^(xi — xq) and yi = 2"^~^^'^Xi, z = 2,3 for transversal coordinates, one has 
obtained |yi| < 2~^ and yi < 2~-^/^, i = 2,3 as for Candes-Donoho-Demanet curvelets. 

Summarizing, the window of a Minkowski curvelets ipA should be a rotated 4d rectangle 
oriented along the four- vector = (1,6^), 

A = n- ({zo2^■/2~n^(l^ 1, 0, 0) + 2-™(0, 0, Z2, zs) + z,2-^/^-^{l, -1,0,0); 

ki-l<Zi<ki + l, i = 0,...,3}), (1.11) 

where the 4- vector (1, 1, 0, 0) is in the direction of propagation, (0, 0, Z2, z^) is orthogonal to 
the propagation plane span((l, 0, 0, 0), ca), whereas (1, —1, 0, 0) is some arbitrary transverse 
direction, here chosen to be the symmetric null direction in the propagation plane. 
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Figure 1: {M = 2). Space-time tiling for a propagation along the first axis. Transversal directions, 
e_L = (0, 0, 1, 0) or (0, 0, 0, 1) are perpendicular to the plane of propagation. The size of the window 
along these two directions is 2~™. 



The following easy adaptation of the construction of the previous paragraph defines a 
frame of Minkowski curvelets satisfying the previous requirements: 



(i) (Fourier scales) start from a one-dimensional Lemarie-Meyer wavelet, ip, and let 

(j > 1), (1.12) 



„ ^/\fT\ ^ f:2/|t2| 



(1.13) 



where ipo is the 'father wavelet' defined in eq. (1.3). 
If ^ S supp (V'™'"') and, say, > 0, then 



\e 



•,2m 



im-j/2 



m « 1^1 « 2-+^v2 



:i.i4) 



as required above. On the other hand, the support of "0™ made up of 'elliptic' 
regions which are irrelevant for geometrical optics where either |^o| ^ 1^1 or |^| <^ |^o|- 
Estimates are very easy in these and we shall not bother any longer about ^™ in the 
remainder of the article. 



(ii) (angular dependence) Let (replacing simply ip^ by ip"^'^ in eq. (1.5)) 



V.-'^-.-(^) = V^-'^-(e)I(|| )i>^>(Map(|| )). 



:i.i5) 



Hence (decomposing ^ G supp (^/;™J'<^) into (.^Oi (e • + ^_l)) there are three Fourier 
scales: the lower Fourier scale, |^o~|^|| ^ ~ ^ ' ^\ ^ 2™"-'/^; the middle one, 
|^_l| ^ 2"^; and the higher one, Co ~ |C| ~ Co + ICI ~ 2'"+j/2. If e = (1,0,0), then 
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Figure 2: (M = 2). Support of a curvelet in Fourier space. Respective scalings are ca • ^ ~ 2™ 
^0 « 1^1 « 2'"+J/2 and < |^| < 2™ for transversal directions. 

supp is included in some enlargement of the rectangle defined by |Ci — Col ^ 

2m-i/2^ 161, 161 < 2"", 161 < 2'"+^/2_ foj. ^j^g real-space window A, the Fourier 
support is oriented along the four- vector eA- It is shifted along that same direction, 
so that ^ supp (■0"*'-''^). 

(iii) The light cone and the removal of the irrelevant zone |6| |€| cut the remaining 
Fourier domain into four connected components, depending on the signs of 6 

161 ~ 1^1- Thus ^'"J'^ may be rewritten as the sum of four functions with disjoint 
supports, which we still denote by ip"^'^'^ to keep notations reasonable. 

(iv) (space localization) Let 

^m,j,e,k^^^ ^ 2-2rn^m,j,e _ (^ko2^ /2-m 1,0,0) + 

2-™(0, 0, k2, ks) + A:i2-J/2-™(l, -1, 0, 0))) . (1.16) 

We let D be the set of all Minkowski curvelet indices (m, j, e, k). If m is fixed, then the 
set B restricts to B"*; if m, j are fixed, then B restricts to B'"'-', and similarly for B"*'^^. 
The Fourier support of A is supp^(A) := supp(V'A) (it does not depend on k). We also 
choose XA € A for each A, for instance 

XA = (A;o2J/2-"^(1, 1, 0, 0) + 2-"^(0, 0, fcs, k^) + fci2-^/2-"»(l, -1, 0, 0)) . (1.17) 

Note that for Minkowski curvelets, 

Vol(supp^(A)) K 2^"^ ^ (Vol(A))-\ (1.18) 
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2 Decay estimates for Minkowski curvelets 



We shall now show estimates for •i/'AiV'A and for the scalar products (V'AiV'A')' showing in 
particular a summable, fast decay in the scaled distance between A and A', in the spirit of 

As in the Candes-Donoho-Demanet construction, we now define Minkowski molecules, 
which have the same real-space and Fourier space decay as Minkowski curvelets. We need a 
good choice of local coordinates for a fixed unitary vector e = e^- Let first (e, e±^i, e±^2) = 
7^^((1, 0, 0), (0, 1, 0), (0, 0, 1)) be an orthonormal basis of M'^, and {x/^,x± = {x±,i, x±^2)) 
the dual coordinates. Then one must choose an adapted basis of vector fields; we let 
V// = Ve = d^o + dx/^; V_L^j = dxj_i, i = 1,2 be a basis of the transversal plane (V_l 
stands for some linear combination of V_L^j, i = 1,2); and Vo// be any derivative in di- 
rections away from the propagation plane, for instance Voff = dx^ — dx^^- Note that 
e - X = x^ — X 1 1 (Minkowski scalar product) gives apparently the 'component' of x along the 
direction (1, — e^) which is 'perpendicular' to the direction of propagation; but both words 
'component' and 'perpendicular' refer to the Euclidean scalar product on M'^, which is to- 
tally misleading in this context. On the other hand, e ■ x is not a (Minkowski) 'projection' 
of X along e either since e • e = 0. In any case, e • x should not be mistaken with e • x, 
which is the Euclidean projection of a; G in the direction e. We also let = ca be 
Minkowski orthogonal to ba, and V*^ = Vg^ 1 1 (^'^ next paragraph, ca will not be 
a null vector any more, and we will have e\ 7^ ba, V*^ h). 

Definition 2.1 (Minkowski molecules) A family (aA)AeD of functions is a family 
of Minkowski molecules of regularity N if 

aAix) = /a (x - 7^'=^ (ko2^/^-"'il, 1, 0, 0) + 2-"(0, 0, k^, k^) + ki2'^l^-'^{l, -1, 0, 0))) , 



(2.1) 

where /a satisfies the following real-space and Fourier space decays, 
(i) Ifr = rii + r_L + rojj < N, 

|(V*//)V/Vl-V:7//(x)|<22-. 

(2™)''^ (2'"+.?'/2)^o// . 

((2™+J'/2eA-x) + (2™|ccx|) + (2™-^'/V°l)) • (2-2) 

The decay is strongest 'along' e^, i.e. for eA • x (away from the direction of propa- 
gation), weaker in the two transverse directions, and still weaker in arbitrarily cho- 



sen directions; in (2.2), \xx\, resp. \x^\ may be replaced by any linear combination 



\x± \ + C|eA ■ x\, resp. + C"|a;_L| + C"\eA ■ x\. 
(ii) (Fourier decay) 

I/a (01 < 2-2™ (((2— +^-/2g^ . ^)) + (2-™!^^!) + ((2-'"-^-/2^o)))"'^ (2.3) 

(recall {{t)) = 1 + t + t~^). Here again, the decay is strongest 'along' ca, weaker 
in the two transverse directions, and still weaker in arbitrarily chosen directions; in 
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(2.3), \$±\, resp. \^^\ may be replaced by any linear combination \ +C\eA resp. 
I^^^l + C"|^_|_| +C"|eA 'CI- Note that |e-^| and |^o| cl^s bounded from below, while 
is not. 

Lemma 2.2 Minkowski curvelets are Minkowski molecules. Furthermore, the following ad- 
ditional Fourier decay holds for centered curv elets. Assume A = (m, j, e, (0, 0, 0, 0)) is 



centered at the origin of space-time. Then (2.3) extends to derivatives ofip^. 



|v7/v1^v:°//^a(OI<2-^™ 



f \ —N 



(2.4) 

Proof. The Fourier decay is a trivial consequence of the prefactor 2"^™ in eq. (1.14) 



and the boundedness of Fourier support of ■0Ai see eq. (1.12) and ( 1.15| ). Fourier tranS' 



forming ip/\ and rescahng variables yields the inverse prefactor 2 ^"^ • Vol^(A) ~ 2^*". The 
differential operator (V*^''// V^-^ V^^y brings down a factor of order (^o - CnY^^ C±Co°^^ 



(2™-j/2)V/(2™)'-i(2™+-'/2)M/ as in eq. (2.2). Finally, if, say, ca = (1,0,0), then 



^i[2'"-j/2(Co-Ci)^^''-2™+j/2Ci(a;i-xO)-2'"(C2x2+C3x3)]^^ ^2 5) 

in rescaled coordinates where all the functions V'(t-), ^(Cq/C^); ^{C)i ^'■'''^'^ (Map((^) have 
bounded, j- and m-independent support in and are C°° with scale-independent || • | loo- 



norms on their derivatives. The real-space decay in eq. (2.2) is now obtained by differenti- 
ating a number of times with respect to ( and integrating by parts. 

Let us finally prove the additional Fourier decay for centered curvelets. Since /ca = 
(0, 0, 0, 0), V'A = 2~'^^ip^'^'^ with m = m^,j = j'A) e = ca. The scaling behaviour of 

W\.iieo/\ei 



is the same as that of ^'"^'-'(0 = V'( 2^" )V'( 2j )• Elementary computations yield 

|V„*"««)| < ma. ( M?^. M!, ii^i'-ilMi' ] < (2.6) 



and similarly 

|Vo//V^"^'^'(OI <max 



22m ' 2i ' 2J i ~ ' ^ 



|VxV'™'^(e)|<max('|ti,^^M^!^)<2- 



22m ' 2i 

More generally, each derivative V, V = ^ ii,V±,Voff produces a supplementary scaling 
factor 2-™-j/2, resp. 2"™, 2-'^+j/2. 

□ 

With this definition, Minkowski molecules (hence in particular Minkowski curvelets) are 
almost orthogonal, in the following sense. 
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Lemma 2.3 (almost-orthogonality of Minkowski molecules) Let (aA)A6D; (a'^)AeD 
be two families of Minkowski molecules of regularity 4N. Then 

\{aA,a^,)\<d{A,A'r'', (2.9) 

where 

d(A, A') := 24|"'l+2|j-y|(i + 1^ _ ^/|2)(i ^ |j - j'|2) 

{1 + dang{A, A') + 4,//(A, A') + d^,±{A, A') + d^^offi^. A')} ; 

(2.10) 

dang{A, A') := 22™™^-''-' )|eA — ga'I is the rescaled angular distance between the Fourier 
supports suppg(A) and suppg(A'), measured in terms of the greater angular span of the 
two; 

dx,//{^, A') := lm+j72<m'+i722™+^/^|eA-(xA-a;A')l+lm'+j'/2<m+i/22"'+^'/^|eA'-(a;A'- 
xa)| is a scaled distance 'along' the direction of propagation; 

'^a;,±(A,A') := lm<m'2'"|eA A {x/^ — ojaOI + lm'<m2'" |eA' A (a^A' — is a scaled 

transversal distance; 

dx,off{A, A') := 2™™('"~-'/^'™' — is a scaled distance in some direction away 

from the propagation plane. 



As in Definition 2.1, dx,± and dx,off may also be redefined by using the coordinates 



along somewhat arbitrary directions. 

Proof. If A, A' G D"^'^ have same scale indices , the n the decay in dang follows from eq. 

(2.3) and the decay in d^j/,dx,±,d^^off from eq. ( [2^ . Note that 2^™ « 1/y^Vo1(A) and 
2^ w 1/y/Vo1(supp5(A)), hence \{aA,a'^,)\ is globally of order 0(1). 

In the general case, the various distances dang,dxj/,dx^±,dx^off should be measured in 
terms of the lower scale, mm{k,k') , where {k,k') = {j /2, j' /2), {m + j /2,m' + j' /2), etc. 
because the window A is not well localized in space (i.e. has rescaled size 3> 1) when 
measured in terms of the width 0{2~^ ) of A if A;' > k. Finally, factors 2~l'"~™ I, resp. 
2-|(m-j/2)-(m'-j'/2)|^ 2-\{ni+j/2)-{m'+j'/2)\ arbitrary order < N may be derived from B, 



yielding in particuliar the required prefactor ^2^1™ ™-'l+2|i _|_ j^T^ _ rn'\'^){l + \j — j 



, N 

i\2\ 



□ 



Lemma 2.4 (decay summability) (i) There exists a constant C such that, for every 
r > 5 and A, A" G D, 

d{A, A')-''d{A', A")-' < Cd{A, A")"^ (2.11) 

A'eD 



(a) (same hypotheses). There exists a constant C > such that, for every A GO, 

Y d{A,A')-'' <C. (2.12) 
a'gd 

Proof. 

We shall only prove (i), since (ii) is simpler and may be proved along the same lines. 
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(i) Assume first that m = m" , j = j", = ba", and sum over A' E ID)"^'J'<^a only. Then 
d(A,A') « 1 + max(2™+^'/2|g^ . {xa - x'^)\,2"'\ef^, ■ {xa-xa')\, 

2'"|ei2 • i^A - XA')\,2"'-^/^\x% - xU) . (2.13) 



Changing the real space localization of A' essentially amounts to shifting one of these 
four scaled quantities by an integer. Hence ^^/gpm,i,e^ (i(A, A')~''d(A', A")"'' < 
Cd(A, A")"*" if and only if the following holds, 

Y,iC + \\x-x'\\r'iC+\\x' -x"\r <iC+\\x-x"\\)-' (2.14) 

x' 

if x, x' G Z^, x' ranges over the integer lattice and || • || is the Euclidean norm on 
M^. This inequality is true if r > 4 for 



[ 

Jl<\\x'-x\\<^\\x"-x 



(1 + ||x - x'||)-"(l + - x"||)-W 



^ (l + ||x-x"||)-^ inf {l + \\x-x'\\ydx' (2.15) 

l<||3:'-a;'||<oo 

and the last integral is convergent. 

(ii) Assume now simply that m = m" and j = j" , and sum over A' E D™'-^. Fix ba, say, 
ba = (1,0,0). Index e^' by a couple of integers I = {hjh), |^iM'2| ^ 2^^"^, so that 
|eA - CA'I ~ 2~-''/2|Z|. Then 

4,//(A, A') + d,,x(A, A') + 4,o//(A, A') « 2'"+^V2|g^ . -xa')\ + 

\l\ (2"^|e^^i • {XA - XA')\ + 2™|e^^2 • (^^A - xa')\ + 2™-J'/2|^o^ _ ^.o^,]^ 

while (ia„g(A, A') ~ 2-^/^|Z|. Note the supplementary scaling factor \l\ in the previous 
expression, a consequence of the 'frustration' due to the non-alignment of ca and 
ca'- However this supplementary factor is not needed for the convergence and may 



be simply forgotten. One reduces to a problem similar to (2.14) but in Z , where the 



two supplementary dimensions account for the angle coordinates li,l2- 

(iii) (general case) Assume e.g. that m < m! and ca = (1,0,0), say. Then the transversal 
region X2 G [(^2 ~ 1)2"™", (/c2 + 1)2"™'] (corresponding to one of the transversal sides 
of A) contains 0(2™ ~™) integer multiples of 2"™ , corresponding to different spatial 
localizations for A'. Summing over all of these cost a factor 0(21™'"*"!). The same 
counting argument goes for the four space-time directions and the angular direction. 
All together one loses a volume factor which is at most 



exactly the same factor found in eq. (2.10). Consider a perpendicular direction 



say, and assume e.g. that m < m" . Summing over all intermediary scales m < m' < 
m" and measuring x±^i is terms of the intermediary scaling 2""^ yields a volume 



15 



factor 21™-"^'! for the sum of (i(A, A') over A'. On the other hand, the scaled distance 
d{A, A") should be measured in terms of the lower scaling 2~"^ and not 2""* , which 



yields the same volume factor in the right-hand side of (2.4). Without the prefactors 
in 1 + |m — m'l^, 1 + \m' — m"p, 1 + |m — m"p in eq. (2.10), one would lose only a 
logarithmic factor 0{\m—m"\ 



to pay for the sum over all intermediary scales. A simple 
exponential factor 2^"l^~^ I would be useless since 2~°l^~^ \ 2-a\x-x \ _ 2-<^\^-^ I and 

the logarithmic factor would remain; but Ylm'=mi^~^\^~''^'\'^)~^i^~^\^' ~^"\'^)~^ ~ 
0(1) provided r > 1. Keeping the assumption m < m" , one must still sum over m' < m 
and m' > m" . If m' < m then (measuring X2 in terms of th e int ermediate scale 2""*) 
a volume factor 0(2!™'"™ I) is lost on the left-hand side of (2.4), and nothing on the 



right-hand side. The scaling factors make up for this loss since 



m'<m 



\m—m'\ 



2—\m—m'\ 



2—\m'—m"\ 



1 + \m — m'p / I (1 + \m' — m 



//\2\ 



< C 



< c 



2—\m—m"\ 

1 + \m — m 
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)\m—m'\ 



m'<m 



2" |m— m"| 

1 + |m — m"P 



(2.16) 



provided 2r — 1 > 0. If on the contrary ml > m" then (measuring X2 in terms of the 



intermediate scale 2 ™ ) a volume factor 0(2 



|m— m'h 



resp. 0(2l™-"^"l) is lost on the 



left-hand side, resp. on the right-hand side of (2.4), all together one has lo st a f actor 
0(2'™' ~™' I) which is compensated by the scaling factors exactly as in eq. (2.16). 



□ 



3 Ultra-violet bounds for Klein-Gordon operator on Minkowski 
space-time 



We consider here the Klein-Gordon operator on flat space-time, 



(/X > 0). General bounds such as the one given in the Main Theorem are greatly improved 
by the use of Minkowski curvelets which are adapted to rather than to the wave operator 
□0. Because has constant coefficients and is so simple, the Klein-Gordon curvelets (iPa) 
we shall be using here are a straightforward modification of those used previously. Up to a 
rescaling, one may assume that = 1; however, in order to keep track of the //-dependence, 
we shall simply assume that /i ~ 1, which allows us to keep the coarse scale at j = 0. One 
then simply substitutes — fj? to in eq. (1.12), so that (1.14) becomes 



\e 



-,2m 



m-j/2 



m+j/2 



(3.1) 



with E{$) :-- 



Lemmas 2.2, 2.3, 2.4 still hold provided one redefines V//, ca and e\ as follows. Let 



be the 'center' of the support of ipA (any point in suppV'A would do equally well). The 



derivative operator V// is defined to be parallel to the hyperboloid ^ 



at ^Ai so 
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that V//(^2 _ ^2 
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0, and contained in the propagation plane. This fixes V// up to 



£ A cA 

normalization. We choose V// := d^^ + ^d^-^ if ca = (1,0,0). Then = (1, |2r,0,0) and 

e\ := (^, 1, 0, 0) is chosen to be Minkowski orthogonal to ca- 

Recall [2] that the retarded Green function of the Klein-Gordon operator is given in 
Fourier coordinates by G'(^) = lime_i.o ^i-f^i-l^sgni^o) ' "^^^ causal Green function is given 
instead by lim£_i.o g2_^2„ig ; and other Green functions are linear combinations of either G 
or with l£(,>o(5(^^ ~M^) and 1£q<o(5(^^ — /^^)- Since Klein-Gordon curvelets are supported 
outside the hyperboloid, the following lemma holds equivalently for any of these Green 
functions. 

We now prove the main theorem, which we restate for the convenience of the reader. Its 
proof follows along the same lines as the previous estimates. 



Theorem 3.1 (ofF-diagonal decay of Green function) Let A G D™'-?', A' G , with 

m,m' G Z, j,j' > 1, m + j /2,m' + j' /2 > 0. Then, for every N >0, there exists a constant 
Cn > such that 



|Ga,A'| < 2-2«"PK™')(i(A,A') 



~N 



(3.2) 



Proof. 

By space-time translation invariance, we may assume that A' = (m', j', e', (0, 0, 0, 0)) 
is centered at the origin of space-time and use the improved Fourier decay (|2.4|). Since 



< 2 2max(m,m') for ^ G supp(^A) nsupp(V'A')) the Fourier part of the decay, i.e. the 

bound |Ga,a'| < 2-2™a^{'^''"')(d£(A, A'))"^, follows as in (i). We shall also prove the space- 
decay in Fourier space this time by using a Fourier transformation, assuming for simplicity 
that the Fourier indices (m,j, ga) and {m',j',e^/) coincide, so that is obtained from 
^A' by a translation in real-space. By definition of G and of the Minkowski curvelets. 



G'a,A' 



11^ 



riV'A'(e)iV(o, 



(3.3) 



where 



(0 = expi2-™ {(A;o2^/2 + ki2-i/^)^^ _ (A;o2^/2 _ ^^2-^/^)^^ _ fc^ . 



(3.4) 



is an oscillation proportional to dx{/S., A'). Then we rewrite Uk{C) as ^-m-'il^l. ' '^^^P- ^-m^^^^ 
(along the two transversal directions) or ^^s+t/II^) and use an integration by parts. Bounds 



for VV'a(Oi ^ = V//,V_L or Vq//, follow from eq. (2.4), yielding precisely the inverse of 
the above scahng factors, 2"™'"-'/^, resp. 2"™, 2~"*"'"-'/^. Similarly 



V 



// 



< 



2-m- 


-i/2 


Vi 






^off 


f 1 







1 

2-m+j/2 



< _ 



< 

~ ^^2 



(3.5) 



17 



For the bound concerning V// we used the following fact (written for ca along the first axis 
of coordinates): 

I'^l/ie - f^') = ^0 - = {Co - C^) + ^(ef - 6) < 2'"-^/^. (3.6) 

These arguments generalize to curvelets A, A' with different Fourier indices by translat- 
ing the curvelet whose window has a larger side. □ 

References 

[1] N. D. Birrell, P. C. W. Davies. Quantum fields in curved space, Cambridge Monographs 
on mathematical physics, Cambridge University Press (1982). 

[2] N. N. Bogolioubov, D. V. Shirkov. Quantum fields, Benjamin Cummings (1982). 

[3] J.-M. Bouclet, N. Tzvetkov. On global Strichartz estimates for non trapping metrics, J. 
Funct. Anal. 254 (6), 1661-1682 (2008). 

[4] E. J. Candes, L. Demanet. The curvelet representation of wave propagators is optimally 
sparse. Comm. Pure Appl. Math. 58 (11), 1472-1528 (2005). 

[5] E. J. Candes, D. L. Donoho. New tight frames of curvelets and optimal representations 
of objects with singularities. Comm. Pure Appl. Math. 57 (2), 219-266 (2004). 

[6] E. J. Candes, L. Demanet. Curvelets and Fourier integral operators, C. R. Acad. Sci. 
Paris, Ser. I 336, 395-398 (2003). 

[7] A. Cordoba, C. Fefferman. Wave packets and Fourier integral operators. Comm. PDE's 
3 (11), 979-1005 (1978). 

[8] Constructive quantum field theory. Proceedings of the 1973 Erice Summer School, ed. 
by G. Velo and A. Wightman, Lecture Notes in Physics 25, Springer (1973). 

[9] C. Fefferman. A note on spherical summation multipliers, Israel J. Math. 15, 44-52 
(1973). 

[10] J. Feldman, J. Magnen, V. Rivasseau, R. Seneor. Construction and Borel summability 
of infrared by a phase space expansion. Comm. Math. Phys. 109, 437-480 (1987). 

[11] J. Feldman, V. Rivasseau , J. Magnen and E. Trubowitz. An infinite volume expansion 
for many Fermions Creen's functions, Helv. Phys. Acta, 65, 679 (1992). 

[12] E. Hernandez, G. Weiss. A first course on wavelets, CRC Press (1996). 

[13] G. Lebeau. Estimations de dispersion pour les ondes dans un domaine strictement con- 
vexe, Journees Equations aux Derivees Partielles, Evian-les-Bains, expose VII (2006). 

[14] V. Mastropietro. Non-perturbative renormalization. World Scientific (2008). 

[15] J. Metcalfe, D. Tataru. Global parametrices and dispersive estimates for variable coef- 
ficient wave equations, jarXiv:0707.1191 



18 



[16] Y. Meyer. Ondelettes, Hermann Editeurs des sciences et des arts, Paris (1990). 

[17] Peskine, Schroder. An introduction to quantum field theory, Addison- Wesley (1995). 

[18] M. A. Pinsky. Introduction to Fourier analysis and wavelets, Graduate studies in math- 
ematics 102, American Mathematical Society (2002). 

[19] V. Rivasseau. From perturbative to constructive renormalization, Princeton Series in 
Physics (1991). 

[20] M. Salmhofer. Renormalization: an introduction. Springer Verlag (1999). 

[21] A. Seeger, C. Sogge, E. Stein. Regularity properties of Fourier integral operators, Annals 
of Math. 134 (1991), 231-251. 

[22] H. F. Smith. A parametrix construction for wave equations with C^'"^ coefficients, Ann. 
Inst. Fourier 48 (3), 797-835 (1998). 

[23] H. F. Smith. A Hardy space for Fourier integral operators, J. Geo. Anal. 8 (4), 629- 
(1998). 

[24] C. Sogge. Lectures on non-linear wave equations. International Press, Boston (2008). 

[25] E. M. Stein. Harmonic analysis. Real-variable methods, orthogonality, and oscillatory 
integral, Princeton mathematical series 43 (1993). 

[26] D. Tataru. Nonlinear wave equations, Proceedings of the International Congress of 



Mathematicians, vol. Ill, 209-220 (Beijing, 2002), available on arXiv:math/0304397 



[27] D. Tataru. Strichartz estimates for operators with nonsmooth coefficients and the non- 
linear wave equation, Amer. J. Math. 122 (2), 349-376 (2000). 

[28] D. Tataru. Strichartz estimates for second order hyperbolic operators with nonsmooth 
coefficients U, Amer. J. Math. 123 (3), 385-423 (2001). 

[29] D. Tataru. Strichartz estimates for second order hyperbolic operators with nonsmooth 
coefficients HI, J. Amer. Math. Soc. 15 (2), 419-442 (2001). 

[30] F. Treves. Introduction to pseudodifferential and Fourier integral operators. Plenum 
Press (1980). 



19 



